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Tuning topological disorder in MgB 2 
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We carried out Raman measurements on neutron-irradiated and Al-doped MgE>2 samples. The 
irradiation-induced topological disorder causes an unexpected appearance of high frequency spectral 
structures, similar to those observed in lightly Al-doped samples. Our results show that disorder- 
induced violations of the selection rules are responsible for the modification of the Raman spectrum 
in both irradiated and Al-doped samples. Theoretical calculations of the phonon density of states 
support this hypothesis, and demonstrate that the high frequency structures arise mostly from 
contributions at q ^ of the E2 9 phonon mode. 
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A few years after the discovery of superconductivity 
in MgPj2, many of its characteristic properties have been 
investigated and mostly clarified. Nowadays a wide con- 
sensus exists about MgB2 beinga phonon-mediated two- 
gap/two-band superconductor 0, 13, 0, 0, 0, where su- 
perconductivity mainly arises from the coupling of the 
holes in the a band with the E2 9 phonon mode, which 
involves the in-plane stretching vibrations of boron atoms 
El El E3 ■ In order to identify the fundamental mecha- 
nisms of the superconducting pairing, a key role has been 
played by a variety of experimental techniques which per- 
mitted the tuning of the electronic and lattice properties 
of MgB?, such as the application of external pressure 
fill Il2j , irradiation-induced disorder Hii. Ilfil , chemi- 
cal substitution (e.g. Al on the Mg-site [lfil ll 7I | and C on 
the B-sitc [HEl). In many of these cases, a reduction 
of the critical temperature, or even a suppression of the 
superconducting phase, have been observed. In doped 
MgB2 compounds, in particular, the debate focuses on 
the relevance and the relative weight of two different mi- 
croscopic mechanisms responsible for the reduction of T c , 
namely the band filling effects which tune the electron 
density of states [H 0, and the disorder-induced 
increase of the inter-band scattering 0, Hj| . 

Since the early investigations on MgB2, Raman spec- 
troscopy has been widely applied because the doubly de- 
generate Eig mode, responsible for superconductivity, is 
the only Raman active among the four optical phonons 
(A 2ll , Bi 9 , Ei„ and E 2g )0,0]- Indeed, the most promi- 
nent feature of the measured Raman spectrum, namely 
the very broad peak (width ~ 300 cm -1 ) centered at 
around 620 cm" 1 [H 0]), has been generally 
identified with the Ei g phonon peak, whose predicted fre- 
quency at the zone center T lies within the 510-660 cm -1 
range H |27| . Despite the almost general consensus on 
the above assignment, however, there remains unsolved 
controversies |2J, |2^, |2Jj . Most of the doubts arises from 
the anomalously large phonon linewidth and from its re- 



markable temperature dependence HHH, IHIlS. Dif- 
ferent mechanisms, based on the strong electron-phonon 
(c-ph) coupling and on anharmonic effe cts, have been 
invoked to explain these features 0, 125L Efl Isil l32l | , al- 
though recent theoretical calculations suggest these ef- 
fects do not fully account for the experimental findings 
|29) . As an alternative idea, the anomalously large line- 
shape was proposed to be related to the excitations of 
phonons outside the T point activated by lattice defects 
and multiphonon scattering |2flj . This hypothesis has 
also been suggested by some exp erimental works where 
the large phonon line-shape [28j, and other weak spec- 
tral features in the Raman spectrum |25j , were ascribed 
to phonon contributions, throughout the Brillouin zone, 
activated by disorder-induced relaxation of the Raman 
selection rules. 

Raman experiments have been also carried out on pure 
and doped MgB2 films [32I I33I l3~3| and on bulk chemi- 
cally substituted MgB 2 compounds [H El EH . Also 
in these cases, the attempts at a quantitative analysis 
lead to ambiguous and sometime contradictory results. 
Contrary to the prediction of one Raman-active mode 
only, the introduction of even small quantities of sub- 
stitutional impurities appears to drive the onset of ad- 
ditional hi gh-f requency (HF) structures around 750 and 
850 cm" 

1 H3, Ea Eli, which are absent or vanishingly 

small in the pure compound. At higher levels of substi- 
tution, the Raman signal further complicates, progres- 
sively approaching the spectrum of the end compound. 
In Mgi_ x Al x B2 a two-steps behaviour is observed: the 
HF structures are well detectable at x = 0.05 and, up to 
x = 0.25, their intensity increases while keeping almost 
constant the frequency positions; for x > 0.25 the inten- 
sities keep growing but the frequencies increase and the 
lineshapes narrow, evolving as a whole towards the spec- 
trum of the parent compound AIB2 fl7l l2*fij. However, 
since C and Al doping induces, at the same time, charge 
doping, volume compression and on-site lattice disorder, 
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it is particularly difficult to discntaglc their effects on the 
Raman spectrum and on the superconducting properties. 

The present paper aims to single out the effect of topo- 
logical disorder by earring out a Raman study of differ- 
ently neutron irradiated MgB2 samples We show 
that neutron irradiation causes the appearance of HF 
spectral structures, similar to those observed in lightly 
Al-doped samples. The comparison of the present spec- 
tra with the literature data allows us to ascribe the unex- 
pected HF Raman features, present in the irradiated and 
lightly Al-doped sample, to a disorder-induced violation 
of the Raman selection rules. Theoretical calculations of 
the phonon dispersion curves support this hypothesis and 
allow us to compare the Raman results with differently 
q- integrated phonon density of states (PDOS). 

Isotopically enriched ( 11 B) polycrystalline MgB2 was 
prepared by direct synthesis from pure elements Three 
identical samples were irradiated by means of thermal 
neutrons at the source SINQ (Paul Sherrer Institut, 
Switzerland) Different irradiation times were ap- 

plied to the samples to accomplish 7.6*10 17 , 1.0*10 19 , 
and 3.9*10 19 n/cm 2 . The most important damage mech- 
anism in our samples is driven by the neutron capture re- 
actions by the residual amount of 10 B present (less than 
0.5%). Lattice defects are created by the recoil of 4 He and 
7 Li, produced by the nuclear reaction, which are emit- 
ted isotropically. This makes the defect distribution very 
homogeneous. After the irradiation, the resistivity pro- 
gressively increases (16, 64, 124 [ifi cm) and T c decreases 
(35.9 K, 24.3 K, and 12.2 K) on increasing the fluencc. 
For sake of comparison low Al-doped samples (5%, 10%, 
20%, 30% with a T c of 36.6 K, 33.4 K, 29.1 K and 24.1 
K) were also prepared following a similar procedure |35j . 

Raman spectra were measured in back-scattering ge- 
ometry, using a confocal micro-Raman spectrometer 
equipped with a charge coupled device detector and notch 
filter to reject the elastically scattered light. The sam- 
ple was excited by the 632.8 nm line of a 30 mW Hc-Ne 
Laser. Raman spectra were collected using a 20x objec- 
tive (laser spot about 10 /xm 2 wide at the sample surface) 
over the frequency range 200-1100 cm -1 , which includes 
the whole phonon spectral region. For each sample Ra- 
man spectra were collected from different points on the 
sample surface and then averaged to avoid possible effects 
of preferred microcrystalline orientation. 

The Raman spectra of the irradiated and Al-doped 
samples are shown in Fig^ Accordingly to literature 
data 0,113, the effect of Al doping on the Raman sig- 
nal is that of driving the onset of a HF spectral struc- 
ture. The Raman spectra of the irradiated samples, also 
shown in FigQ] present an unexpected strong dependence 
on the fluence, here observed for the first time. Quite 
interestingly, the two series are surprisingly similar: in 
both cases, HF structures appear and become dominant 
at highest values of Al doping or neutron fluence. The 
same features, occurring at about the same frequencies, 
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FIG. 1: Background subtracted Raman spectra of Al-doped 
(a) and neutron irradiated (b) MgB2 samples. Spectra are 
shifted upwards of a constant value. 



have been observed also in light C-doped MgB2 |32j 
and disordered films |33j ]. Since neutron irradiation has 
the only effect of disordering the lattice, these results 
suggest that the HF structures in the chemically substi- 
tuted samples are produced mainly by the lattice disorder 
as well. In this context, the weak hardening of the fre- 
quencies of the HF structures, well detectable only for 
30% Al-doped samples, could be likely ascribed to the 
Al-driven lattice compression. 

In light of these results, and because of the unclear 
origin of the HF structure, the core questions arise: what 
is the origin of the HF structure and what is the role of 
the topological disorder to make them visible? And, most 
importantly, is the q = E2 ff phonon mode the one really 
probed in the Al-doped and irradiated samples? As we 
are going to discuss, we believe that the answer to the 
last question is negative and that we are actually probing 
some kind of phonon density of states. 

In Fig. [2]the Raman spectrum of the highest irradiated 
sample is compared with the PDOS of MgB2 measured by 
inelastic neutron scattering Since no selection rules 
apply to neutron scattering, the excellent agreement be- 
tween the experimental data over a wide frequency range 
(500-1000 cm -1 ) suggests that in disordered MgB2, the 
Raman selection rules are notably relaxed and the Ra- 
man spectrum basically reflects the PDOS. On the other 
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FIG. 2: Raman spectrum of the most irradiated sample (full 
line) compared with the PDOS from the neutron experiment 
of Ref. pjj (squares) and with the calculated E2 g projected 
PDOS (dashed line). All the reported quantities are empiri- 
cally normalized to the same peak intensity. 
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FIG. 3: Left panel: phonon dispersion of MgB2 within the FC 
model. The thickness of the lines in the phonon dispersion 
reflects the magnitude of the Ei a component for each mode. 
Right panel: the corresponding total PDOS (dashed line), 
and the E^g projected PDOS (shaded area). 



side, the apparent disagreement below 500 cm suggests 
a remarkably different weight of the phonon modes in the 
Raman and the neutron experiment. In the present case 
we can assume that the topological disorder, induced by 
neutron irradiation or by chemical substitution, mostly 
tunes the violation of the Raman momentum selection 
rules but not of the eigenvector selection. 

For a more quantitative check of this empirical hypoth- 
esis, we calculated the phonon dispersion of MgB2 within 
a force-constant (FC) shell model. Four elastic springs, 
specified by their tensor connecting B-Mg, Mg-Mg first 
neighbors, and up to B-B second neighbors were con- 
sidered. The corresponding twelve parameters were ob- 
tained by fitting the first-principle calculated phonon dis- 
persion from Ref. j2||. Since the FC model is meant to 
reproduce the bare phonon dispersion, in the fitting pro- 
cedure we did not include the E 2g phonon branch along 
the r — A direction, which is known to be strongly af- 
fected by the el-ph interaction. The effects of the el- 
ph interaction involving the E 2g boron mode with the 
electronic a band were included through the self-energy 
renormalization of the phonon frequencies, by means of 
the two-dimensional Lindhardt function, which, for the 
parabolic almost two-dimensional a bands, results to be 
a very good approximation. Specific details, as well as 
the fitting values of the twelve FC parameters, can be 
found in Ref. ■ 

In Fig. |3 we show the so-calculated phonon disper- 
sion curves for the different phonon branches /i and 
the total PDOS (dashed line). Our calculations enable 
to evaluate the E 2g component of each phonon mode, 
f/q,/f = |£q,;j • ££ 2a | 2 ; which we believe to be effectively 
probed in our Raman measurements, by simply project- 
ing the generic phonon eigenvector iq^ on the E 2g one, 
6e 2 ■ Thus, in Fig. ||we show also the E 2g charac- 
ter of each phonon branch, as quantified by the thick- 



ness of the lines, as well as the -E2g-projected phonon 
densiy of states (i?2 9 -PDOS), shown as shaded area in 
the right panel. The resulting i?2 9 -PDOS, smeared by a 
Lorentzian function with width T = 100 cm" 1 , to sim- 
ulate the finite phonon lifetime, is also shown in Fig. [21 
for comparison purposes. The good agreement between 
the i^g-PDOS and the Raman spectrum of the highest 
irradiated sample supports the idea that the topological 
disorder induces a violation of the Raman momentum se- 
lection rules, so that the effective spectra correspond to 
q-integrated quantities. The slight difference on the low 
frequency side could be ascribed to secondary contribu- 
tions, disorder-activated, from other Raman non-active 
modes, such as the B\ g . It is also worth to notice that 
the peak in the PDOS at about 580-600 cm" 1 and the 
HF structures at ~ 800 cm" 1 , shown in Fig. [3J arise 
from the contribution of different regions of the Brillouin 
zone, respectively from the q-states close to the T-A and 
to K-M. 

This observation suggests that the onset of the HF 
structures in the disordered samples could reflect the tun- 
ing of the breakdown of the momentum selection rules as 
induced by the amount of disorder. In this perspective 
we model the Raman intensity by 

where N c is the number of sampling points in the Bril- 
louin zone and where the momentum cut-off q c depends 
on the amount of disorder and describes the q-space ef- 
fectively probed by the Raman measurements. In partic- 
ular, q c — > for ordered and pure compounds, thus se- 
lecting only the Raman q = signal, whereas for highly 
disordered samples q c is larger than the size of the Bril- 
louin zone and the total q-integrated PDOS is probed. 
On a physical ground, we expect that also the phonon 
scattering rate, T qc , is highly sensitive to the amount of 
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FIG. 4: Raman intensity from Eqs. Q for different val- 
ues of the cut-off q c : (from the bottom to the top) q c = 
0.30,0,45,0.60,0.75 a- 1 . Curves are shifted upwards for clar- 
ity. 



confirmed by the direct comparison between the Raman 
spectrum of the most irradiated sample and the neutron 
PDOS We also showed that the loss of the Raman 
momentum selection rules, tuned by the amount of topo- 
logical disorder, can explain in a natural way the onset 
of the HF structures. Thus the present results provide a 
clear and unambiguous explanation of the Raman spec- 
tra in MgB 2 , and address Raman spectroscopy as a very 
sensitive technique to investigate the role of topological 
disorder in MgB2 and, more in general, in systems where 
a strong coupling between the lattice and the electronic 
degrees of freedom exists. These findings also shed a new 
light on the interpretation of the previously reported Ra- 
man spectra in MgB2 based compounds. 



disorder, here parametrized by q c , with r ge smallest for 
q c — > and increasing with increasing q c . To take into 
account these effects we assumed a linear dependence on 
q c , namely T qc = aq c , where a = 120 A/cm -1 . 

In Fig. 21 we show the intensities calculated using Eq. 
□ for different values of q c (0.30, 0,45, 0.60, 0.75 A" 1 ). 
For |q| < 0.3 A -1 we are essentially probing only the fiat 
dispersion region close to the T point shown in Fig. |3J 
and only the low energy structure at ~ 600 cm' 1 is vis- 
ible. However, for larger q c the loss of the momentum 
selection rules is strong enough to cause probing also of 
the |q| states close to the M and K points (we remind 
that |q M | = 0.39 A -1 and |q K | = 0.68 A" 1 ), making the 
HF structure of the S 2g -PDOS of Fig. rapidly visible. 
Finally, for q c larger than the size of the Brillouin zone, 
all the |q| states are equally probed, and the signal I(oj) 
approaches the full E^-PDOS (Fig. [2J|. The direct ex- 
perimental vs. theoretical comparison in Fig. [21 together 
with the remarkably similar trend shown by the spectra 
in Fig. Hand the calculated intensities in Fig. 0] reveal 
that, also in disordered systems, the Raman spectrum 
mostly originates from the Ei g phonon branch. Giving 
the key role of the Ei g mode, this finding explains the 
frequently observed correlation between the Raman sig- 
nal and the sup erconducting properties of MgB2 based 
compounds [H El El El Ej 

In conclusion, in this Letter we presented Raman spec- 
tra from differently neutron irradiated and Al-dopcd 
MgB 2 samples. We have shown that, on increasing the 
neutron irradiation dose (i.e. topological disorder) as 
well as Al-doping (i.e. topological disorder plus charge- 
doping), HF structures appear and rapidly become the 
dominant spectral features. The remarkable similarity 
among our spectra and those reported in literature from 
lightly C-doped samples suggests that the topological dis- 
order, which is unavoidably present in chemically substi- 
tuted as well as in irradiated sample, is mostly respon- 
sible for the modifications of the Raman signal. This is 
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